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Abstract 

We compute the Andreev reflection coefficient of a flux of thermal quasiparticles in ^He-B incident 
upon various two-dimensional vortex conflgurations. We flnd that, for the same number of positive 
and negative vortex points, the reflection coefficient is much reduced if the points are arranged in 
pairs (which corresponds to a gas of vortex rings in three-dimensions) rather than random (which 
corresponds to a vortex tangle in three-dimensions). This results is consistent with measurements 
performed by Bradley et al. at the University of Lancaster. 



PACS numbers: 

67. 30. em Excitations in He3 

67. 30. he Vortices in Hc3 



Experimental and theoretical studies of turbulence in superfluid ^He and ^He have re- 
vealed features which are similar to what happens in the classical turbulence but also features 
which have no classical counterpart. The study of quantum turbulence at very low temper- 
atures is particularly interesting because, in the absence of the normal fluid, the turbulence 
consists simply of a tangle of reconnecting, thin vortex filaments. Since all filaments have 
the same core radius and the same circulation, the problem of turbulence is reduced to the 
geometry and the topology of the tangle. Unfortunately, unlike what happens in the study 
of ordinary turbulence, only few methods of visualization are possible in liquid helium. In 
superfluid ^He-B, the Andreev scattering technique can be used to detect vortex filaments. 
This technique, pioneered at Lancaster, is based on the fact that the energy dispersion curve, 
E = E{p) of quasiparticles of momentum p is tied to the reference frame of the superfluid, 
so, in a superfluid moving with velocity v^, the dispersion curve becomes E{p) + p ■ (see 
the review article Thus a side of a vortex line presents a potential barrier to oncoming 
quasiparticles, which can be reflected back almost exactly becoming quasiholes; the other 
side of the vortex lets the quasiparticles to go through. Quasiholes are reflected or trans- 
mitted in the opposite way. The vortex thus casts a symmetric shadow for quasiparticles at 
one side and for quasiholes at the other side j^, and, by measuring the flux of excitations, 
one detects the presence of the vortex. 

However the extrapolation from the scattering of quasiparticles off one vortex to the 
scattering off many vortices may be nontrivial. In our earlier work we used a simplified 
2-dimensional vortex point model, and found that the Andreev shadow caused by simple 
configurations of several point vortices is not necessarily equal to the sum of shadows of 
individual, isolated vortices (in the cited work such a phenomenon was called a 'partial 
screening'). This result may have non-trivial implications for the interpretation of Andreev 
reflection measurements of the vortex line density in turbulent '^He-B. 

The present work is motivated by the experimental observations of Bradley et al. ^ of 
the transition from a gas of vortex rings to a dense vortex tangle. Still using the simplified 
2-dimensional model, the question which we ask is whether a gas of vortex rings and a dense 
tangle have a different Andreev signature. 

We start with describing the motion of vortices which, in the considered two-dimensional 
approximation, become vortex points. Such a two-dimensional system of point vortices in 
the inviscid fluid is known as the Onsager's point vortex gas. Each vortex point moves 
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with the flow fleld generated by all other vortices. In the (x, ?/)-plane, the i^^ vortex point 
located at time t at the position rj(t) = Xi{t)i + where i and j are, respectively, the 

unit vectors in the x- and ^/-directions, generates the fluid velocity fleld 



where Ki is the circulation generated by the i vortex. In the considered two-dimensional 
point vortex gas modelling turbulent ^He-B, = ±k with k = Trh/m ^ 0.662 x 10^^ cm^/s, 
where k is the quantum of circulation in ^He-B, m is the mass of ^He atom, and the signs 
plus and minus denote respectively a vortex generating the anticlockwise rotation of the 
fluid in the {x, ?/)-plane, and an antivortex generating the clockwise rotation. At the point 
r, the velocity fleld created by the system of vortices is given by the superposition of the 
velocity flelds generated by all vortices. 

An isolated pair of a vortex and an antivortex can be considered as a two-dimensional 
model of a three-dimensional vortex ring. Such a pair moves through the fluid, in the 
direction orthogonal to the line connecting the vortex and the antivortex, with the velocity 
V = K/{2Trd), where d is the distance between the vortex points in this pair. Obviously, 
in the point vortex gas a vortex and an antivortex can be considered a pair (a 2D vortex 
ring) if d is much smaller than the average distance between the center of this pair and the 
nearest vortex which does not belong to this pair. 

To model turbulence in ^He-B, we consider the following three configurations of the vortex 
points: 

1°. A spatially random system of A^ point vortices of the same polarity; such a system has 
net circulation = Nn. This configuration can be considered as the two-dimensional 
model of a polarized vortex bundle 

2°. A random, statistically uniform system of A^/2 vortices and A^/2 antivortices in the case 
where the distance between any two vortex points is not much smaller than L/-\/iV, where L 
is the size of the computational domain, so that vortices and antivortices are not organized 
in pairs. In such a system there is no net circulation, Kf = 0. This configuration can be 
considered as a two-dimensional model of a vortex tangle. 

3°. A system of vortex-antivortex pairs such that the distance d between the vortex points in 
each pair is the same, but locations and orientations of pairs are random. This configuration 
can be considered the two-dimensional model of a gas of vortex rings. 




(1) 
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The configurations 2° and 3° can be characterized by two geometric quantities, d and the 
average distance, say a between the centers of two neighbouring vortex-antivortex pairs. The 
relevant non-dimensional parameter characterizing the configuration is ( = d/a < 1. Gas of 
small, compared to the mean intervortex distance, vortex-antivortex pairs (configuration 3°) 
corresponds to C ^ 1- Of particular interest is the transition between configurations 3° and 
2°. Such a transition corresponds, in a modelling sense, to the experimental observations of 
Bradley et al. ^\ of the transition from the gas of vortex rings to the dense vortex tangle. 
In the considered two-dimensional system of vortex points this transition can be modelled 
by increasing the parameter C, from small values corresponding to d of the order of lOO^o, 
where ~ 0.85 x 10^^ cm is the zero-temperature coherence length, to values of the order 
of unity. 

We consider the two-dimensional problem of ballistic propagation of quasiparticle excita- 
tions in the flow field of the point vortex gas and the gas of vortex-antivortex pairs described 
by configurations 1°, 2°, and 3°. Neglecting spatial variations of the order parameter, in the 
presence of the flow field the energy of thermal excitation can be written as 



E=^ep2 + A2 + p.v,(r, t), (2) 

where ep = p^/(2m*) — ep is the kinetic energy of a thermal excitation of momentum p 
relative to the Fermi energy ep ~ 2.27 x 10^^^ erg (here and below the numerical values 
are taken at zero bar pressure [0]), m* ^ 3.01m = 1.51 x 10^^'^ g, with m being the mass 
of the ^He atom, and Aq = l.TGhsTc ~ 2.43 x 10^^^ erg, with ks being the Boltzmann's 
constant and Tc the critical temperature, is the superfiuid energy gap. (It should be noted 
that the superfiuid energy gap was assumed a constant value, Aq. because we are concerned 
with the behaviour of thermal excitation at distances r from the vortex core larger than the 
zero-temperature coherence length ,^o-) Excitations with > are called quasiparticles, 
and excitations with < are called quasiholes. 

Following the approach of Refs. 3, 3] we assume that the interaction term p • varies on 
a spatial scale which is larger than so that the excitation can be regarded as a compact 
object of momentum p = p(t), position r = r(t), and energy E = E{p, r, t) given by Eq. ([2]). 



Using the method developed in Refs. 



Eq. ([2]) can be considered as an effective, semi- 



classical Hamiltonian yielding the following equations of motion: 



^ dEjp r, t) ^ ep ^ 

dp ^ep^ + A^m* ' ^' ^ ^ 



4 



P--^M._||p., (4) 



Introducing the non-dimensional variables [sj 



R= ^, V, = — v„ r = -, n = — , = — , (5) 



where R = (X, Y), pp = ^Jlm*tp ~ 8.28 x 10 g cm/s is the Fermi momentum, and 
to = ^oPf/Ao ~ 2.9 X lO^S, in the non-dimensional form the Hamiltonian ([2]) and the 
equations of motion ([3])-(llD become 

E = VA2(n2 - 1)2 + 1 + m*m-\^U ■ V,(R, r), (6) 

and 

R= , ^ ^ n + — vrV, R, r , (7 

n = -m*m-V2V(n- V,(R, r)), (8) 

where now A = dA/dr, and we introduced the non-dimensional parameter X = ep/ ^o- In 
the following calculations we assume A = 10^. 

We solve numerically the equations ([7D-([H]) of ballistic motion of quasiparticles in order 
to calculate the propagation of thermal flux, generated by the point source of thermal 
excitations, through the point vortex gas configurations 1°, 2°, and 3° introduced above. 
The numerical method is described in our earlier work 3|]. 

In turbulence experiments in '^He-B, the properties of the vortex tangle or the gas of vortex 
rings can be studied by measuring the heat which is transported by thermal excitations 
through the velocity field of the vortices 9| (see also the review article Ref. l|). A net flux 
of excitations (and, hence, energy) results in the case where there is a (small) temperature 
gradient. In this case the heat carried by excitations generated by the source (and, therefore, 
incident on the vortex gas) is 

SQ^nc= I NpVFE^^^dTdE, (9) 

where 5T <^ T is a temperature difference between the source of thermal excitations and 
the opposite side of the system, Np = my/2m*eF/ (vr^/j,'^) is the density of states at the 
Fermi energy, vp = Pf/^* ~ 5.48 x 10^ cm/s is the Fermi velocity, and f{E) is the Fermi 
distribution, which, at the ultra low temperatures, becomes the Boltzmann distribution 
f{E)=eM-E/kBT). 



In the case of an isolated vortex, one side of it presents a potential barrier to oncoming 
quasiparticles, which are Andreev reflected as quasiholes almost exactly back to the source. 
In our previous work ^ we found that the Andreev shadow, i.e. the maximum distance from 
the vortex core past which a quasiparticle with the kinetic energy ep is not Andreev reflected, 
is 5*0 = 37r(Ao/ep)^ in our dimensionless units. In the subsequent work [3] it was shown that, 
due to partial screening, the Andreev shadow of simple configurations of several vortices is 
not necessarily equal to the sum of shadows of isolated vortices. Hence we have reasons 
to expect that in a gas of point vortices or in a gas of vortex-antivortex pairs (modelling, 
respectively, a vortex tangle or a gas of vortex rings) partial screening will strongly affect 
the heat flux carried back to the source by Andreev reflected quasiparticles. 

The numerical simulation of the heat flux reflected by the point vortex gas was carried out 
in the rectangular domain. At t = the vortices are randomly distributed within the square 
sub domain whose non-dimensional coordinates are —5 x 10'^ < X < 5 x 10^,-5 x 10^ < 
y < 5 X 10^ (the size of this domain corresponds to the experimental estimates 0,1^, l3, 11|). 



The heat source is located at (—1 x 10^, 0) so that the angle, if between the X-axis and 
the beam of quasiparticles varies between — 7r/4 and ir/A. The flux generated by the source 
is modelled by = 25000 quasiparticles whose non-dimensional initial energies, Hq {1 < 
Hq < 1.688) and directions, ipo are uniformly distributed. Since the source of excitations 
is located sufficiently far from the vortices, the initial momentum of quasiparticle whose 
initial energy is Hq was calculated from Eq. (E]) as Hq = [X^^Hq - 1)^^^ + 1]^^^- For the 
Boltzmann distribution, f{E), introducing E = E/ksT the integrand in ([9]) reduces to 
the r(3) distribution which, in our numerical calculations, was generated by the standard 
subroutine and the resulting values were discarded if < Aq/ ksT. In the typical low 
temperature experiments, Ao/(fcBT) ^ 10; this value was used throughout all calculations. 

A trajectory of each quasiparticle was found by numerical solution of Eqs. ©-([H]) for 
random initial energies (momenta) and directions of motion. Having identified the trajec- 
tories of quasiparticles that are Andreev refiected (as quasiholes) by the vortex gas or a gas 
of vortex-antivortex pairs, for a particular realization, R of the initial configuration of the 
vortex gas the reflection coefficient is then calculated as 

^ K K 
ff = -y^aiH,o with qo = y^H,o, (10) 

U 

where aj = 1 if the j**^ quasiparticle is reflected, otherwise aj = 0, and K = 2.5 x lO'*. 



FIG. 1: (Color online) Zoomed view of configurations 1° (left), 2° (center), and 3° (right); + (red) 
- vortices, Q (black) - antivortices. 




FIG. 2: Reflection coefficient as a function of the total number of vortex points. Lines, from top 
to bottom, correspond to configurations 1°, 2°, and the gas of vortex-antivortex pairs 3° with 
d = 100^0- 

This procedure has been carried out for A^r reahzations of the initial configuration of the 
vortex gas, and the ensemble average reflection coefficient was calculated as fr = (/^) for 
configurations 1°, 2°, and 3°. Since the error of ensemble averaging decreases with the 

— 1/2 

number of realizations as A^j^ , to achieve a reasonable (few per cent) accuracy we used up 
to N^i = 2.5 X 10^ realizations. 

At this point it seems useful to show a zoomed view of each of these configurations, 
see Fig. [1] Fig. [2] shows the reflection coefficient, fr as a function of the total number of 
vortices, A^ in the system. As can be seen from this figure, in all cases the refiection coefficient 
increases with the total density of the vortex points. The reflection coefficients in the system 
1° of vortices of the same polarity and in the vortex gas 2° with zero net circulation are 
close, although in the latter case the partial screening seems to have a slightly bigger effect. 




FIG. 3: Reflection coefficient, for 13 vortex-antivortex pairs (A^ = 26), as a function C = d/a, 
where d is the distance between the vortex and antivortex in a pair, and a is the average distance 
between the centers of nearest pairs. 



However, in the case where the vortices of opposite polarities form pairs ('rings'), partial 
screening plays a far more pronounced role and the reflection coefficient falls by almost an 
order of magnitude (in the calculation illustrated by the bottom line of Fig.[2]the dimensional 
distance between the vortex points in each pair was assumed d = lOO^o)- 

As was already mentioned above, of particular interest is a behaviour of the reflection 
coefficient during the transition between configurations 2° and 3°. Such a behaviour can be 
characterized by the reflection coefficient, as a function of the non-dimensional parameter 
(. This function is illustrated by Fig. [3] for = 26 (for other values of the behaviour of fr 
with ( remains qualitatively the same). As can be seen, in the considered example, with 13 
vortex-antivortex pairs, the increase of ( = d/a from ( = 0.01 (gas of small vortex-antivortex 
pairs) to ( = 0.2 (almost random, disordered point vortex gas) is accompanied by almost 
flve-fold increase of the reflection coefficient, from fr ~ 3.5 x 10^^ to 0.158, respectively, the 
latter value being not very different from that for the random point vortex gas (/,. ~ 0.169). 
Therefore, the results shown in Fig. [3] seem to indicate that the transition, observed by 
Bradley et al. j4|, from the gas of vortex rings to the relatively dense vortex tangle can 
be detected by the significant (nearly an order of magnitude) increase of the coefficient of 
reflection of the heat flux carried by thermal excitations. 

The experiment 4] was performed at temperature T = 0.16 and pressure P ^ bar 
where Tc is the critical temperature. Quantized vortices were generated by an oscillating 
grid and detected by two vibrating wires placed near the grid. A beam of quasiparticles 
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illuminated the grid. In the presence of vortices a fraction of the quasiparticles is Andreev 
reflected, reducing the damping of the vibrating wire; this damping is caused by the asym- 
metry of the quasiparticles and quasiholes incident upon the wire. The transient response 
of the fractional change of the damping was measured as a function of the velocity of the 
oscillating grid. It was found that at high grid velocity (4.5 to 7.8 mm/s) the fractional 
reduction of damping was from 0.25 to 0.4 and recovered slowly in 10 to 15 s, whereas at 
small grid velocity (1.9 to 2.9 mm/s) the fractional reduction of damping ranged from 0.025 
to 0.1 and recovered quickly in less than 1 s. The experimenters suggested an interpretation 
based on the recovery time, i.e. that at high grid velocity quantum turbulence is created 
which slowly decay and disperse away. On the contrary, at small grid velocity the vorticity 
is in the form of a gas of vortex rings no larger than 5 /xm, which quickly move away (the 
translational velocity of a vortex ring is inversely proportional to its size). Our computed 
results are consistent with this interpretation: we found that the reflection coefficient of a 
gas of vortex- antivortex points (the two-dimensional equivalent of a gas of vortex rings) is 
much less than that of a random vortex gas (the two-dimensional equivalent of a vortex 
tangle), which agrees with the experimental observation. 

We arc grateful to D. I. Bradley, G. R. Pickett, S. N. Fisher, W. F. Vinen, and also to the 
participants of the Workshop on Topics in Quantum Turbulence (Abdus Salam International 
Centre for Theoretical Physics, Trieste, March 2009) for stimulating discussions. 
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